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We study the effect of site diagonal disorder on the pairing amplitude 
by a perturbation method. Using an extended Hubbard model with the 
intersite attraction we analyze fluctuations of order parameter in presence 
of non-magnetic disorder and discuss the instability of various solutions 
with p-wave paring. 
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1. Introduction 

In case of an anisotropic superconductor non-magnetic disorder leads 
eventually to pair breaking effect via Abrikosov-Gorkov formula [1, 2, 3, 4, 
5], like for magnetic impurities, for any value of a coherence length £. By 
analogy to an isotropic s-wave pairing, for superconductors with an enough 
large coherence length £, the amplitude of paring potential |A(r/)| tends 
be the same for all bonds but disorder induces the effect of fluctuations in 
the amplitude leading, in this way, to destruction of pairing. Here we will 
calculate fluctuating potentials standard deviations < <5|Aj.,'| 2 > and < ef > 
Their ratio T =< <J|Ay| 2 > / < ef > will be a criterion of pairing potential 
fluctuations [6]. 

2. P— wave solutions for a clean system 

We start a single band, extended, Hubbard model with effective nearest 
neighbour intersite attraction W%j (Wu = 0) [7]. Taking the Fourier trans- 
form of the Green function for a clean system we can write its equation of 
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Fig. 1. (a) The Fermi surface and (b) the Free energy 5F as a function of pair- 
ing potential |A| (|A| = \J A^. + A^) for different solutions: dipole (D) real (R) 
and complex (C) order parameters, respectively; for 7 band of Sr 2 RuC>4 electronic 
structure (t' = 0.45i, n/2 = 0.66). The intersite attraction W=-0.446t, tempera- 
ture T = 



motion: 
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where A^ = A x sin k x + A^ sin k y defines 2x2 matrix of a pairing poten- 
tial. For a clean system tij can be expressed in fc-space: e fc = Y^j t a e ~ xRiik 
= — 2t(cos fe^+cos ky)—At' cos k x cos k y , where t represents the nearest neigh- 
bour site amplitude of electron hopping, while t' corresponds to next near- 
est neighbour one. For the assumed solution [8] A& = io y o ■ d(k), with 
d(k) = (0, 0,cP(k)) and d z (k) = The corresponding free energy F for 
a finite temperature T can be found from the standard formula: 

k 

where E k denotes the eigenvalue. To perform numerical calculations we 
have fitted our one band system parameters to the realistic 7 band struc- 
ture of Sr2Ru04 [5, 9, 10, 11] Fig. la presents the corresponding Fermi 
surface. For the above assumptions we have found three solutions with 
p-wave pairing. Namely, depending on relative values A x and A y , the 
dipole one (A x ^ and A y = 0), the real one (A x = A y ) and the com- 
plex one (A x = i A y ). They correspond to minima of free energy curves 
(SF(\A\) = F(\A\) -F(0)) in Fig. lb, denoted by D, R and C, respectively. 
The interaction parameter used in calculations (W=-0.446t) were chosen to 
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give T c = 1.5K as for clean Sr2Ru04. One can easily see that the complex 
solution reaches the global minimum of free energy F. 



3. Fluctuations of pairing potential. 

In this section we investigate the stability of superconducting p-wave 
states in presence of a weak disorder. Here, we apply the same strategy as 
in [6] and we treat random site energies £j as perturbations. To proceed 
we write the Dyson equation for a Green function G(i,j;u>) evaluated at a 
frequency uj: 

G(i,j;co) = G ( i ,j;^)+^G°( i ,j»F i G(/,j», (3) 

i 

where Vi = Ei<T3 is the impurity potential matrix. 

Following Eq. (3) we express quantity order parameter Ay, in the low- 
est order of perturbations by means of disordered Green function and 
calculate the mean square deviation of the paring Ajj along the bond of 
nearest neighbour sites i and j. For site independent energies £»: 

<<5|A ij | 2 >=r ij <e?>. (4) 
Finally, we calculate the coefficient r y - [6]: 
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Having found the pairing potentials (Fig. lb, Eq. 1) for the clean 
system, we calculated V (Eq. 5) for all three solutions we obtained. In 
all cases V has a very small value (of order ~ 10 -8 ). This implies that 
fluctuations of Ay are relatively small in this system. Interestingly, a real 
type solution is characterized by the smallest fluctuations (r = 5.15 • 10~ 9 ) 
while T = 13.93 • 1(T 9 for the dipole solution and T = 8.12 • 10~ 9 for the 
complex one. This could mean that the real solution is favoured by disorder. 



4. Conclusions and Discussion 

We have analyzed the effect of a weak disorder on a p-wave supercon- 
ductor in context of newly discovered superconductor Sr2Ru04 [11, 12]. 
Unfortunately the order parameter structure in this compound is still un- 
known [12]. Fitting our system parameters to its 7 band structure we have 
asked about the stability of various solution in presence of disorder. We 
have found three solutions with the same critical temperature Tc = 1.5K: 
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dipole and real solutions with line nodes in the gap and complex one with 
a finite gap in any direction. Note that, all these solutions have the same 
critical temperature Tc and the Abrikosov-Gorkov formula [2, 4] does not 
differentiate any of them. Our preliminary results at T = indicate that the 
complex type of solution has the global minimum of free energy F but the 
real one is favoured by disorder. That result was obtained in one band model 
in the lowest order of perturbation method and should be confirmed by a 
more sophisticated method like Coherent Potential Approximation [4, 5] 
considering more realistic three bands structure of S^RuC^. 
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